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Abstract 

The setups for precise measurements of the time structure of Nuclotron internal and 
slowly extracted beams are described in both hardware and software aspects. The CA- 
MAC hardware is based on the use of the standard CAMAC modules developed and 
manufactured at JINR. The data acquisition system software is implemented using the 
ngdp framework under the Unix-like operating system (OS) FreeBSD to allow the easy 
network distribution of the online data. It is demonstrated that the described setups are 
suitable for the continuous beam quality monitoring during the experiments performed 
at Nuclotron. 
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Introduction 


The experimental program at Nuclotron for heavy ion collisions mm, few body [5] and 
polarization HIS] physics requires good quality of the internal and extracted beams. The 
time structure of the beam plays a crucial role in the experiments requiring high inter¬ 
action rate. One of the parameters reflecting the beam time structure is the coefficient 
itidc defined as 
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t2 — ti 


-I 2 


{dN/dt)dt 





(dN/dt)'^ dt 
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where dN/dt is the beam current, {t 2 — ti) = Text is the time of the beam exposition. 
The real experiment requires that the coefficient should be about 0.85—0.95 for the 
efficient data taking. 

Another problem of the experiment related with the beam time structure is the correct 
estimation of the data acquisition dead-time. This is important not only for the cross 
section evaluation, but also for the polarization measurements due to effect of the dead¬ 
time distortion [^. For instance, large dead-time of the data acquisition (DAQ) due to 
oscillating interaction rate could significantly reduce the asymmetry values during the 
measurements of the analyzing powers [Hi- Such effect is especially important for the 
polarization measurements requiring the absolute normalization of the beam intensity or 
interaction rate in the wide range as in the TPD experiment on the measurements of the 
induced tensor polarization of deuteron beam travelling through matter Emo]. 

The Nuclotron internal target station (ITS) is an efficient tool to perform nuclear 
physics studies at energies from hundreds of MeV to several GeV per nucleon with both 
polarized and unpolarized beams m- However, in this case the useful events yield is the 
time dependent function of the beam and moving internal target interaction. Recently, 
new system of the ITS control and data acquisition has been developed m- In particular, 
the internal target mechanical vibrations were studied in details using this new system 

m- 

The Nuclotron extracted deuteron beam time structure has been studied previously 
during data taking of the TPD experiment Ellin], on the measurements of the soft 
photons yield [13] , and by the BM@N experiment [T] . The CAMAC and VME [TI] based 
DAQ systems were used for the measurements performed at 4V El na na and 6v [n 
beamlines, respectively. These measurements were performed in the trigger mode of the 
DAQ systems, when the information was read for each trigger. However, such mode 
cannot provide correct information on the beam time structure at high trigger rate due 
to dead-time distortion. These circumstances motivate us to organize the setup with 
negligible dead-time and to use it for the Nuclotron beam time structure measurements. 

The goal of the paper is to report the details of the setups for the Nuclotron internal 
and extracted beam time structure measurements as well as first results obtained with 
deuteron, ^Li, and beams. The paper is organized in the following way. The chapterjl] 
contains the details of the hardware and software of the setups for the Nuclotron beam 
time structure investigations. The results on the extracted beam time structure are 
presented in chapter |2j The results on the investigation of the internal beam-target 
interaction are reported in chapter jS] The conclusions are drawn in the last chapter. 
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Figure 1: CAMAC hardware. See text for description. 


1 Setup hardware and software 

1.1 CAMAC electronics 

Each setup contains input register, two identical scalers, and pulse generator CAMAC 
modules. Their functional scheme is outlined in Fig.[I] All CAMAC units were developed 
and produced at JINR. Both setups have the same modules disposition in the read-out 
crate. 

The input register collects three inputs — accelerator cycle begin and two triggers 
— which lead to the LAM assertion. The input register is the only LAM source in the 
read-out crate, and the KK009 crate controller m is programmed to set the interrupt 
request (IRQ) at each LAM occurrence. Each of this IRQ leads to the invocation of 
the CAMAC interrupt handler function from the spill loadable kernel module (see 
chapter Ol) . Each of two triggers corresponds to scaler module should be read-out at 
current interrupt handler invocation. The two triggers should alternate strictly, so the 
Spill DAQ is able to handle unproper trigger bit combinations and trigger sequences. 

The KV005 [12] generator was used to produce the trigger signals sequence with the 
base frequency of 10"'^ Hz, which corresponds to the time slice of 100 /rs. 

The beam-target interaction intensity is measured by monitor scintillation counter, 
whose signal is split to fed the 0th input of both scalers. The odd scaler counts during 
odd time slice and is read-out during even one, and the even scaler — vice versa. This 
scheme allows us to avoid undercounts, only up to 1 pulse could be erroneously dropped 
or added in each time slice. 

Because we read from scalers the least 16 bits only at each 100 ^s time slice, the 
setups are able to consume counts with rate up to ~ 6.5 x 10® Hz. This is more than 
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Figure 2: Overall Spill DAQ layout. See text for description. 


enough because scalers are able to input up to ~ (0.7 1) x 10® Hz. 

1.2 Spill DAQ system 

Both setups have very similar software however each uses own hardware, in particular 
the main host. The overall Spill DAQ layout is pictured in Fig. [21 where we can see the 
main host under the FreeBSD OS as the rectangle entitled “Spill DAQ”, and optional 
hosts — as the rectangles “visualization client”. This main host deals with the CAMAC 
hardware using the PK009 ISA adapter for the KK009 CAMAC crate controller [15] . 
In the OS kernel context it has the camac(4) facility [17], spill(4) CAMAC kernel 
module, as well as the modules in the netgraph(4) package [H] style (Fig. |3|). Fig. |3| 
shows the data streams transportation from the spill(4) CAMAC kernel module through 
the nodes, which are represented as rectangles with node name (up), type (left), and ID 
(right). The data streams are propagated along the graph edges, which are formed by 
pairs of connected hooks [18] (octagons in Fig. [S]). The data are injected into graph by 
the cEunacsrc node and demultiplexed by the fifos node into more than one identical 
streams, which could be transferred both locally and remotely. In Fig. |2]we can see two 
of the data stream consumers: writer( 1) and b2h( 1), whose standard inputs are fed by 
standard outputs of ngget(l) instances, which reads the fifos local outputs through 
the netgraph sockets provided by ng_socket(4) node instances. Each fifos output as 
well as their input could be connected and disconnected without disturbing other one(s). 
Due to fifos loading during the OS boot sequence all load user commands (see Tabled]) 
are independent from each other. The “visualization client” in Fig. [2] executes one of the 
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Figure 3: The Spill DAQ core is implemented by the ngdp graph. See text for description. 


possible clients of the b2h(l) server — the ROOT script clients. C . 

The netgraph(4) package [TB] and ngdp framework [1^1 US] involvement reduces 
the DAQ system implementation efforts essentially, because a number of software mod¬ 
ules are ready to use: ng_camacsrc(4), ng_fifos(4), ngget(l) [20], writer(l) [21] . 
ng_ksocket(4), ng_socket(4). The produced data could be distributed in form of 
both the ngdp packet stream(s) (for raw cycle-by-cycle ones) and ROOT framework 
PB] TMessages (for b2h( 1) produced histograms). Specifically for the Spill DAQ we im¬ 
plement only the spill(4), spillconf(l), spilloper(l), and b2h(l) software modules. 

The spill(4) module is intended to work with the corresponding CAMAC hardware 
(see fTT]) . It complies with requirements of the camac(4) and ng_camacsrc(4), so im¬ 
plements the CAMAC interrupt handler function. This handler recognizes the following 
interrupt occurrences (events): 

• begin of burst (BoB), 

• odd and even trigger at end of 100 ps time slice. 

In total two data packet types are produced to contain as follows: CYC_BEG — the BoB 
timestamp, CYC_END — the read-out scaler data as the sequence of 50000 16-bit little- 
endian values to cover up to 5 s of burst. The CYC_END packet type is produced at 
obtaining the 50000th end-of-slice interrupt. 

The CAMAC hardware description and handling are separated from the spill(4) 
module’s source and grouped together in the single spilLhardware.h header file. In the 
current implementation it uses macro interface kk(9) specific for the KK009 crate con¬ 
troller m instead of the crate controller independent one camac(9), because the former 
interface allows us to perform each CAMAC cycle slightly faster [53]. 

The spill(4) module can be conhgured at startup by the spillconf(8) utility and 
controlled during operation by the spilloper(8) one. To operate with the Spill DAQ 
it has the configuration file (named by default $NGDPHOME/etc/spillsv.conf) in the 
Makefile format (see also make(l)). This file establishes the correspondence between 
the user commands (“targets” in make(l) terminology) and actions which should be 
performed. This hie is textual and could be revised easily. So typical user command look 
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Table 1: The user commands to control over Spill DAQ. 


Command 

Description 

load 

loads and configures the spill(4) and ng_ccmiacsrc kernel mod¬ 
ules, connects latter to ng_fifos input 

unload 

counterpart for previous 

loadw 

loads the writer( 1) module and connects it to ng_f if os output. 
RUNFILE and DATADIR variables could be defined 

unloadw 

counterpart for previous 

loadb2h 

loads the b2h(l) module and connects it to ng_fifos output. 
RUNFILE variable could be defined 

unloadb2h 

counterpart for previous 

continue 

starts the handling of CAMAC interrupts 

swstart 

starts the handling of CAMAC interrupts with software BoB im¬ 
itation by callout(9) mechanism 

pause 

counterpart for two previous 

init 

initializes the CAMAC hardware 

finish 

counterpart for previous 

clecuiall 

resets all histograms in the b2h(l) 

saveall 

stores all b2h(l) histograms into ROOT TFile on hard disk (if 
RUNFILE variable was supplied for loadb2h) 

disconn 

disconnects all clients from b2h( 1) 

status 

outputs status messages from present modules through sys- 
logd(8) 

seelog 

starts the syslogd(8) output file qdpblog viewing by tail(l) 


like the following: 

make -f spillsv.conf RUNFILE=test loadw 

All the defined user commands we summarize in the Table [TJ 

The b2h(l) (for “binary to histograms”) module is intended for ROOT histograms 
filling from the packet (binary) data representation obtained from ng_f if os node output. 
Due to simple nature of our setups we use the compiled-in histograms conhguration from 
hconf.h header file. The 50000-bin histograms for at least last cycle and sum during 
current run, and optionally their Fourier transforms are booked at b2h(l) startup. At 
each CYC_END packet arrival these histograms and some statistics are updated. Namely, 
for 50 bins of 100 ms the mean of counts, its standard deviation, dispersion and dispersion 


Table 2: The time structure characteristics for the Nuclotron slowly extracted beam. 
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Figure 4: The monitor counts in the 0.1 s window (left) and Fourier transform of the 
full burst (right) during the single cycle from the slowly extracted deuteron beam of 
Tkin = 3.5 A-GeV in F6 focus. See text for description. 


Table 3: Nuclotron internal beam-target interaction characteristics. The time limits for 
the Aide (see equation O) calculations are: = 1000 ms, ^2 = 3500 ms. The frequency 
ranges for the mechanical vibrations of the internal target (Af,Ts) are 12..14 Hz for C 
and 13..15 Hz for GHz. 


Beam, 

(A-GeV) 

Target 



AN{Af) /o/'i 

O.SAftot 

for A/ as / 1../2 (Hz) 




0..1 
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5.4 
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0.35 

d, 3.5 

C 

0.918 30.1 

0.9 


5.9 

1.1 

0.53 


per mean are calculated and printed to standard error output. The b2h(l) listen(2)s 
on port 12342 for the client registration. For each registered client the b2h(l) sends 
all histograms by TMessages once per cycle. Also all configured histograms are written 
optionally to ROOT TFile periodically or at b2h(l) termination as well as at SIGHUP 
signal obtaining. The SIGINT signal resets all histograms, SIGUSRl and SIGUSR2 switches 
to read-and-discard mode and vice versa. For user termination in accuracy manner the 
SIGTERM signal should be used. 

The ROOT script clientS.C comvects to b2h(l) server, obtains THIF histograms en¬ 
coded into TMessages, and draws each of them in separate TCanvases. To allow the full 
user interaction with these histograms the script should be terminated (i.e., by twice 
<Ctrl><C> pressing), because the proper coexistence of two endless loops (for TSocket 
reading and XI1 [53] events handling) could be managed only in the executable like 
histGUI(l) [20]. 
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Figure 5: The monitor counts in the 0.1 s window (left) and Fourier transform of the 
full burst (right) during the single cycle from the slowly extracted ^Li beam of = 
3.5 A-GeV in F4 focus. See text for description. 

2 Nuclotron extracted beam results 

The first setup version was implemented in 2009 and successfully used during 7 Nuclotron 
runs in 2011-2015 to measure the extracted deuteron, carbon, and ^Li beams time struc¬ 
ture. The scintillation counter (100 x 40 x 4 mm^) with XP2020Q PMT movable by 
motor was situated in the beam halo near the F4 focus of the VPl transport beamline. 
The PMT pulses shaped by 4F-115 discriminator were used as the input monitor counts 
for the setup. 

The measurements were performed during the data taking on the measurements of 
the soft photons yield m- The typical averaged number of the counts in the 100 fis time 
slice was 28, 24 and 7 per beam spill for the ^Li, deuteron and carbon beam, respectively. 
In Figs. 0]—[5] the typical monitor counts proportional to the beam-target interaction rate 
during a single cycle are shown for the beam of deuterons, ^Li, and carbon with 3.5 A-GeV 
kinetic energy. The abscissa is a time bin number, the ordinate is the number of counts 
in the corresponding bin. The bin width is equal to the time slice, 100 fxs. The typical 
0.1 s time windows are shown on the right panels instead of the full burst timescale. The 
Aide value (see equation ([1])) could be obtained easily from these data in online or offline. 
For the single cycles depicted in Figs. 01—[5] the Aide during the full burst time — G) 
are represented in Table [2l We are not show the whole run (per-bin cumulative sum) 
histograms because the single cycle statistics is large enough. The Fourier transforms 
of the corresponding full burst statistics are shown in Figs. 0]— [5] on the left panels. As 
one can see, some frequencies (50 Hz and harmonics) are emphasized. Their origin is a 
time structure of the slowly extracted beam. The ratio values for some set of 

frequency ranges A/ = /2 — /i (in Hz) are represented in Table 0] The Aide coefficient 
is smaller than 0.5. This reflects the large width of the frequency spectrum (see right 
panels in Figs. 0HB]). On the other hand, the contribution of the harmonics proportional 
to 50 Hz is lower than 7 % in the worst case of ^Li. One can conclude, that significant 
improvement of the extracted beam time structure is required for the high intensity 
experiments, for instance, by the feedback from the monitor counter to the accelerator. 
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Figure 6: The monitor counts in the 0.1 s window (left) and Fourier transform of the 
full burst (right) during the single cycle from the slowly extracted beam of T^in = 
3.5 A-GeV in F6 focus. See text for description. 


3 Nuclotron internal beam results 

The second setup version for the Nuclotron internal beam was implemented in 2014 
and successfully used during 2 Nuclotron runs (June’2014 and March’2015) to measure 
the time scans of the beam-target interactions in the accelerator ring. The scintillation 
counter (100 x 100 x 4 mm^) with FEU-30 PMT observes the beam-target interaction 
point of the Nuclotron internal target station (ITS) [T2] from 1 m distance under 60° angle 
in horizontal plane to produce the input monitor counts for both ITS and current setup 
DAQ systems. The analog signal from PMT with the rising time of ~ 5 ns and total 
duration of ^ 15 ns is shaped by the 4F-115 discriminator with ^ 40 mV threshold. The 
duration of the output pulse is ^ 100 ns. In the second setup version the counting is 
started by hardware at the same time offset for each cycle (with Ri 50 ns uncertainty) 
after the BoB event occurrence. This allows us to work with the whole run histograms 
because the subtle counting start variation in time can not “smudge” the time structure 
on such additive picture. 

The quality of the electronic chain has been checked using ^°®Ru radioactive (RA) 
source. The Fourier transform of the monitor counts from the ^°®Ru RA source is shown 
on the left panel in Fig. [7] for the whole run. These data were obtained while the 
150 A-MeV deuteron beam was accelerated and circulated in the Nuclotron ring, however 
the internal target was not run. The averaged number of counts in the 100 /iS time slice 
is ^ 15. The data demonstrate white noise picture that proofs the quality of the used 
electronic chain. The right panel in Fig. [7] shows the Fourier transform of the monitor 
counts from both the deuterons of the same energy on CH 2 target and the ^°®Ru RA 
source simultaneously. The frequencies corresponding to target mechanical vibrations 
and to the n ■ 50 Hz harmonics are clearly seen. 

The data on the deuteron beam time structure were obtained simultaneously with the 
data taking for the DSS experiment [3] with 10 /rm CH 2 foil and 8 /rm carbon multiwire 
targets with the typical intensity of 5 8 • 10® deuterons per spill. The typical averaged 

number of the counts in the 100 /iS time slice was ~ 15 and ~ 180 per spill at the beam 
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Figure 7: The Fourier transform of the monitor counts from the ^*^®Ru radioactive source 
without (left) and with (right) Nuclotron internal deuteron beam of = 150 A-MeV 
on the CH 2 target for the whole run. See text for description. 


energy of 250 A-MeV and 3.5 A-GeV, respectively. 

In Figs. [5] and IHl we can see the full burst timescale (left panels) of the monitor 
counts during the whole run (per-bin cumulative sum) from the beam of deuterons with 
250 A-MeV kinetic energy interacting with the CH 2 and C targets, respectively. The 
Fourier transforms of the corresponding full burst statistics for the CH 2 and C target 
are shown in Figs. |8] and [9] on right panels. The same for the 3.5 A-GeV deuteron 
beam is shown in Figs. ITOl and fTTl One can see some frequencies (50 Hz and harmonics) 
other than originated from the mechanical target vibrations are enhanced, however their 
relative intensity is decreased with the frequency increasing. Their origin can be a space- 
time structure of the Nuclotron internal beam. The ratio values for some set of 

U.oJVtot , 

enhanced frequencies / in narrow ranges A/ = / ± 1 Hz are represented in Table |3l The 
“A/its” Column contains this ratio for mechanical vibrations frequency of the internal 
target (~ 13 Hz for C and ss 14 Hz for GH 2 in the represented data). The results confirm 
the earlier observation [12] that the mechanical vibrations are larger for GH 2 target. The 
iFdc coefficient is ^ 0.9 at 3.5 A-GeV for both targets. The contribution of the n • 50 Hz 
harmonics is only about 8 10 %. However we observe the signihcant difference in ATjc 
coefficient and contribution of the n • 50 Hz harmonics for CH 2 and C targets at low 
energies. For carbon target the ^ 0.8 and the n • 50 Hz harmonics contribution is 
only 1.5 %. Therefore, one can conclude that the time structure of the internal beam 
interaction with carbon target allows one to perform experiments at ITS. For CH 2 target 
the Aide is only 0.56 and the n • 50 Hz harmonics contribution is ^ 15 %. The reason of 
this difference is the subject of further studies. 


Conclusions 

The results of the present work can be summarized as follows: 

• Two similar setups for precise measurements of the Nuclotron beam time structure 
are designed using standard GAMAG modules produced at LHEP and DLNP. They are 
commissioned for the experiments at internal target and at the slowly extracted beam. 
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Figure 8: The monitor counts during the whole run (left) and their Fourier transform 
(right) for Nuclotron internal deuteron beam of r^in = 250 A-MeV on the CH 2 target. 
See text for description. 


Both setups are working with the 100 p,s time slices, so they are able to investigate the 
frequency domain up to Ri 4.5 kHz. 

• The data on the beam time structure could be easily distributed online through network 
on the cycle-by-cycle basis due to Spill DAQ design on base of the ngdp framework [191 [20]. 
Therefore, these setups are suitable to work continuously during the whole Nuclotron run 
and provide data for the beam quality estimations. 

• The first setup version was successfully used to measure the time structure of the slowly 
extracted beam of deuterons, ^Li, and carbon at 3.5 A-GeV. 

• The second setup version was used to measure the time structure of the beam-target 
interaction inside the Nuclotron ring at ITS mm- It has been found the influence of 
the mechanical vibrations of the internal target on the data taking. 

• The developed setups can be efficiently used as a tool for the permanent monitoring of 
the Nuclotron beam quality during the experiments performed at Nuclotron. 
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Figure 11: The monitor counts during the whole run (left) and their Fourier transform 
(right) for Nuclotron internal deuteron beam of = 3.5 A-GeV on the C target. See 
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